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<D ■ Abstract 
Q 

' The quasi-channeling of positive and negative relativistic particles in a bent crystal is studied 

using the classical deflection function. It was shown that the potential scattering in a central field 

Qh| of bounded ring-like potentials may produce the "reflected" and "refracted" fractions of scattered 

Q ■ particles. For particles with positive charge the "reflected" fraction is mainly presented; at the same 

. time we predict that for particles with negative charge the "refracted" fraction should dominate. 

O , 

* ^ I The effect of "empty core" for central scattering is also discussed. The average deflection angles for 

>> 

r~| ' volume "reflection" and "refraction" are derived for accepted potential model of the crystal. The 



calculated average "reflection" angle is in satisfactory agreement with recent experimental data 

> 

QQ ■ PACS numbers: 61.85.+p; 29.27.-a; 45.50Tn; 41.85-p; 25.75.Gz 
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INTRODUCTION 



In the early 80-s, the experiments [2, y, |J] have discovered the effect of volume capture 
of relativistic particles in a channeling motion by uniformly bent crystal. Several attempts 
have been made to find a specific mechanism of capture (see e.g. [sl) because the coherent 
collisions in a central field can not cause the channeling. Taratin and Vorobiev studied the 
passage of high energy particles in a bent crystal using the computer simulations and 
found that the particles (both positive and negative) are deflected to the convex side the 
bend of the atomic planes by an angle of about 2 * 9c {6^ - critical angle). They called 
this a volume "reflection" by bent atomic planes. The significant beam divergence and 
multiple scattering of the quasi-channeled particles in the crystals prevented observing this 
effect for almost two decades. The experiment jl, 8] of 70GeV protons scattered by a short 
bent crystal is the first explicit demonstration of the volume "reflection". It should also 



br the 2001 run 



ll| by the same 



be mentioned that preliminary results of crystal collimation at RHIC 
demonstrated some interesting results which were recently explained Id . 
effect. 

My purpose here is to work out the analytic formulas for volume "reflection" and related 
effects, and discuss the principle difference between deflections of positive and negative 
particles. The standard continuous screening potential for the bent crystal is assumed to be 
a good approximation. The small angle formula for classical scattering 12| are not applicable 
here, because the continuous potential could not be considered weak near the turning points 
although the angles of deflection are small. This article will deal with exact solution of 
potential scattering by specific model of the crystal's potential. 



MODELS OF POTENTIAL FOR BENT QUASI-CHANNELING AND CHANNEL- 
ING 

The charged particles incident on a bent single crystal along the tangent directions to 
the crystallographic planes, experience a collective field from atoms of the tangent areas of 
these planes. So, the starting point for studying the quasi- channeling and channeling in 
bent crystals is the interaction potential for binary ion-atom collisions. 
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Uatomir) = ^/(^), (1) 

where Zp^^ are the bare nuclear charge of the projectile and crystal atoms, Rt-f is Thomas- 
Fermi (or Firsov) screening length and / is a screening function, 

i 

The screening function approaches 1 in the limit r and decays exponentially for large 
r. This is only possible if = 1- Different analytic forms of the screening function 

/ are in use. The Moliere approximation!^, /m, uses {a^} = {0.1,0.55,0.35}, = 
{6.0,1.2,0.3}. Other choices include the Ziegler-Biersack-Littmark potential[l^ fzBL with 
{a,} = {0.1818,0.5099,0.2802,0.0281}, and {A} = {3.2,0.9423,0.4028,0.2016}. Now, if 
we find the average potential over a cylinder with the radius R and the root-mean-square 
temperature displacement u of atoms in gaussian form, we receive the following expression 
for continuous cylindrical potential 

Pi '^t^T-F tiT-F y IRt-F 

1^ l^+exp(-A^^^5 ^)er/c(A-^- ^^^^))- (3) 



V2u Rt-f V2Rt-f V^u 

Here p is polar coordinate from center of cylinder, U = ttZpZtRt density of 

atoms in the bent plane, erfcQ is the complementary error function[15]. 

For charged projectiles , the screen potential of atoms ([1]) over all space is positive or 
negative depending on the sign of Zp. It follows that the average potential of each plane Up{p) 
and the whole cylindrical crystal potential U{p) = ^ Up{p) will also be positive or negative 
depending on the charge of the projectile (Fig.l (A)). In the intervals < p < R — Nd, 
R < p < oo the potential U{p) decreases to zero. In the interval R — Nd < p < R, potential 
does not have a singularity and is smooth almost periodic function with some constant 
positive or negative pedestal. The value of pedestal is quite significant. It is 13.3 eV for 
Si < 110 > plane (T=300K) with maximum peak of the potential 35.9 eV, and 10 eV for 
Si < 111 > with maximum peak 34.7 eV, i.e. about 1/3 of maximal value of the screening 
potential. 

When the particles enter the bent crystal from outer side, the influence of the pedestal 
might be significant because the outer potential ring has bigger magnitude. Most deflection 



experiments are done using a different geometry when the particles enter the piece 

of crystal from the front end and the influence of the pedestal in this case is quite small. The 
value of pedestal can be subtracted and 'effective' potential and its rectangular model are 
shown in Figs. 1(C), (D). However, the signs of scattering potentials still remain the same. 
To estimate the trajectories of particles and deflection function we will use the following 
parameters for the crystal Si < 111 > used in[l|: Uo = 25eV, a = 0.78A, d = 3.13A (see 
Fig.l(D)). 

CLASSICAL DEFLECTION FUNCTION ON RING POTENTIAL 

The path of a particle in a central classical field is, as well known, symmetrical about a 
line from the center to the nearest point in the orbit (OA for positive particles in Fig.2(A)). 
Hence the two asymptotes to the orbit make equal angles ip with these lines. The angle 
X = 11 — 2ip through which the particles are deflected as they pass the center of the potential 
varies in the interval — tt < x < tt when the angle ip varies in the interval < (/? < tt. For 
relativistic particles, the angles (p are near 7r/2 and x ~ 0, so it is convenient to introduce the 
angle a = 7r/2 — (see Fig. 2 (A), (B), (C)) which is positive when the particle deflected to 
the left in Fig. 2(A) and negative when the particle deflected to the right. At high energies. 



there is no "orbiting" or "spiral" scattering 



16| where the angle if can have a larger interval. 



However, there is always a possibility of different signs for deflection angle x = 2a in the 
case of a single ring or periodic system of rings concentrated near a circle, as we will discuss 
below. 

We shall first consider the deflection of a particle of mass m and energy E passing the 
boundaries of a positive rectangular ring Uo > with the outer radius R and thickness a 
(see Fig. 2 (B)). The path of a particle is refracted by the angle (6*0 — 6*1) and (/5i — Pq) 
to the incident directions, at the point of entrance and exit of the ring potential. These 
angles define the deflection angle a = x/2 for over-barrier (quasi-channeling) motion of the 
particles: 

a{b,E) = {9o-9,)-{(3,-(3o), (4) 

with the nearest point A on the path located inside the ring. The tangent continuity of a 
particle's momentum gives the relativistic relations between pair of angles do, 61 and /5o, Pi- 
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^1 - (l))sin^{9i) = sin\9o) - 0, (5) 
= sm2(/5o)-0, (6) 



where the 

2EUn 



jt;2 _ ^ig. 



2„4 



(7) 



is the squared Lindhard's angle = ±9^ with the sign of the potential. In deriving 
(IS]), (El), we made one assumption, 2EUo » Uq, which is, in general, satisfied for nonsin- 
gular potentials and wide range of energies including the nonrelativistic. From geometry of 
scattering (Fig. 2 (B)) and relations ([5]), ([5]), we can receive the explicit expressions for angles: 



52 / 52 

9q = arcsin{\l 1 - — ), 9i = arcsinU 1 - _ ^^^J (8) 



/3o = arcs^ni,|l- j^^), A = arcs^n{^l - ^^^^^). (9) 

as functions of impact parameter 6, energy of particle and parameters of potential i?, a, 
Uq. From now on, we will use the normalized parameters: 

b = b/R, a = a/R, d = d/R, ba = b/{R - a) = b/{l - d). (10) 

The defiection angle (jl]) then can be written: 



p. . bWi-b^-^i-^-b^) . K{^i-bl-^i-<l>-bl) 

aib, E) = arcstni , — arcsmi , , 11 

where the first term correspond to (^^o — ^i); and second to {(3i — (3q) in Eq.(j4]). 

It is worth to note that the Eq. ( ITTl) was derived for the particles transmitted inside the 
ring, but it is, in general, valid for other cases as well. In the following, we consider the 
modification of this expression. 



Case 2P. If the impact parameter b satisfy the inequality ^/T^^7i){l — a) < b < (1 — a), 
the last radical in second term of (ITTl) vanishes. The defiection function (fTTj) in this interval 
will be 

a{b, E) = arcsin[ p=^= ) — arcsin[ -= — ). (121 

V 1 - V 1 - 
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It monotonically increases from negative values to positive (Fig. 3 (A)). 

Case 3P. If the impact parameter (1 — a) < 6 < y/1 — (f), the particle does not pass 
across the inner boundary of ring and the whole second term in (|TT1) vanishes. The deflection 
function (ITTi) in this case continues to increase (Fig. 3 (A)), 



b{\/l-b^ - yl-(l)-b^) 

a{b,E) = arcsin{ ), (13) 

VI - 

becomes the well known formula for the scattering by a spherical potential 'well' or cylindrical 
disk. 

Case AP. For a narrow interval of impact parameters ^Jl — (j) <b<loTR — b<^, 
the particles are reflected from the outer boundary of potential and this effect is responsible 
for the volume "reflection" . The second radical in (11 II) vanishes in this case and deflection 
function is given by 



^Ca/i - b^) 

a{b, E) = arcszn(-^^^==^). (14) 
VI - 

In this same manner we can study the trajectories of particles scattered by a negative 
rectangular ring Uo < (Fig. 2 (C)), and it gives the exactly same expression ffTTj) for 
deflection function. The only difference, that (p is negative, produces just three cases in 
comparison with positive ring. 



Case 2N. For negative potential VI ~ > 1 the impact parameter in the interval 



(1 — d)<6<(l — a)Vl ~ 0; the first radical of second term in (fTTj) becomes zero. The 
modified deflection function. 



(h F^ bWi-b^-^i-<p-b-) . b^{^i-<p-bl) 

Vl - Vl - 

monotonically decreases from positive values to negative (Fig. 3 (C)). This interval of impact 
parameters corresponds to the reflection from inner wall of the negative potential. 



Case 3N. If the impact parameter lies in the interval (1 — a)Vl — < 6 < 1, the whole 
second term of (ITTl) becomes zero. The particle freely moves inside the negative ring and the 
form function in this case has the expression ( fT3i) of Case 3P . It was pointed out already 



that this gives the formula for classical scattering by a spherical potential well 
deflection function is shown in Fig.3(C) by the doted line. 



121. This 
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The sequence of considered cases can be violated if the radius of the ring R > = . 
Indeed, in the Case 3P, where y/1 — (f) < 1, the interval (1 — a) < 6 < ^/l — (f) vanishes 
if R increases and d for positive potentials. In the Case 3N where a/1 — (p > 1 the 
interval (1 — d)y/l — (p <h<\ vanishes if R increases and a ^ for negative potentials. 
In both cases, the interval of transmission the particle without the the touching the inner 
wall disappearres, and the deflection functions have discontinuances (Fig.3(B),(C)). This 
situation takes place in experiments Q], where E = 70GeV, (j) = 0.58* 10~^ (Si < 111 >), 
2a = 1.56 * 10^^°m, and Rc = 28cm. In periodic structure, which we consider below, 
the potential thickness, a, should be substituted by the period of structure, d. This gives 
Rc = = Im which is also less than radius of the bend R = 1.7m. 

EFFECT OF "EMPTY CORE" 

There is interested feature of the classic scattering in central fields which we would like 
to discuss here. If the particles are scattering by a positive potential disk (i.e. a field with 
U = for p > R and U = Uo ioi p < R) , the deflection functions (see Eq. (fT3l) ) are strictly 
positive (doted line in Fig.3(A)) and monotonic in the interval < 6 < ^/l — (p, i.e. the 
"disk" scatters strictly as a repulsive potential in this region of impact parameters. We have 
the opposite situation for scattering by negative potential "well" (i.e. a field with U = 
for p > R and U = —Uo for p < R), where all scattering angles are negative (doted line in 
Fig.3(C)). The "well" scatters strictly as a attractive potential. Now we cut off the central 
circular core of these potentials with radius R — a to receive positive or negative "rings" . In 
the positive central ring, all the particles with impact parameters < b < a/1 — 0(1 — a) 
are deflected with a negative sign (Fig.3(A)), i.e. empty core for positive ring scatters as 
an attractive potential for widest range of impact parameters. In the negative central ring, 
all the particles with impact parameters < b < (1 — a) are deflected with positive sign 
(Fig.3(A)), i.e. empty core here scatters as a repulsive potential. 

The physics of this effect is simple. The inner wall of the positive ring creates a force 
toward the center of the ring and particles obtain the negative scattering angle. The inner 
wall of the negative ring creates a force in upward direction and the scattering angle for 
penetrated particles obtain the positive sign. However, to some extent the behavior of the 
particles in the area of "empty core" is independent of the potential form, i.e. different 
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forms of the potential can produce the same effect. We can observe that a projectile with 
some impact parameter b and fixed point of exit (point E in Fig.2(B,C)) has already defined 
direction of motion, no matter what shape of the potential constitutes the ring. Indeed, 
the momentum p and angular momentum pb of projectile must conserve. These imply that 
the path of the projectile must come out from point of exit E and be tangent to the circle 
of radius b in the area of "empty core" (Fig. 2(B),(C)). Hence, the point of exit and the 
impact parameter enables us to obtain the scattering angle of particle without knowing the 
the exact shape of potential inside the ring. There is another useful feature of this effect. It 
produces the coherent action in wide interval of impact parameters greater than the period 
of ring, so it might be useful for applications. 



REFLECTION AND REFRACTION BY SYSTEM OF RINGS 

Finally lets consider the scattering by the periodic system of the potential rings mentioned 
in Sec. The relativistic particle moves with impart parameter b consequently 'refracted' on 
entering and leaving each potential ring. It is easy to understand, that the deflection function 
for this case will be the sum of the deflection functions of the individual rings (fTTj) : 

N 

a{b,E)=J2c^^Hb,E). (16) 

i 

Here 



a^'''{b,E) = arcsin{- 



. bU^l-bl-^^l-<l>-bl) 

arcstn[ j, (17) 

VI - 



where the normalized impact parameters are 



I — ta 1 — a — ta 



In the relativistic case all pairs (W 1 — 6^ — y 1 — — 6^), (y 1 — — y 1 — — 6| 



are 



extremely small and substitution arcsin{x) x can be done in (fTT]) - (fT71) without losing the 
accuracy, e.g. the Eq. (fT6|) can be rewritten 
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"^^'""^^^^ Tr^l Tr^l • ^ 

If the impact parameter b<{l—Nd), all terms in Eq. fll6l) . (HM are present and the particles 
transmit into the empty core of potential (quasi-channeling). The deflection angles in this 
case were discussed in Sec. If the impact parameter is such that the following condition 
occurs 

b> (l-kd) (20) 

for some specific index k, all terms of f|T9l) with indexes i > k vanish and only external rings 
with indexes i < k contribute to the scattering angle. Overall, Fig. 4 shows the deflection 
angles for positive and negative particles. The deflection functions are continuous functions 
of impact parameter for a certain interval of energies E > or a certain interval of ring's 
radius R < j^, which illustrated in Fig. 4(A), (C). The opposite case, 

R> Rc=^, (21) 



which we discussed in Sec. is also plotted in Fig. 4(B), (D). In this case, the deflection 
functions are discontinuous, i.e. the incident coherent beam is divided by crystal's period into 
two separate parts. We can call them "refracted" and "reflected" parts, because "refracted" 
particles are mostly deflected to the convex side of ring and "reflected" particles deflected to 
the opposite side of the bend, although it does not correspond to what is called "refracted" or 
"reflected" in optics. These parts are separated by a gap depended on potential shape. For 
our ring's model, the gap can be estimated from Eq. ffTOj) . "Reflected" and "refracted" parts 
have a distribution of angles rather than deflection at a single angle. In fact, the formula 
for distribution ( fT9l) allows to estimate the average angle for "reflected" and "refracted" 



fractions. For positive particles, this is an average integral over the interval {l — id)^/l — (p < 



b < {1 — id)y/l — (j). Because the deflection function is almost periodic, we can take second 
non distorted interval (Fig. 4(B)) with i = 1: 



(l-a)^W 6(Vl - 62 - ./l - - 62) 



{d — a) VI — J{i-d)^/r^ yl 
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>^ (22) 

VI - 

where 6a, 6i are defined in ( ITOi) . This integral can be estimated, but the result is too cumber- 
some to reproduce it here. However, if we disregard the small terms ,(P jCp"^ ,dd,(pd,(pd,etc. , 
the final expression for average reflection angle a+ will be: 

a+ = -J^ (203/2 _ 272ci3/2 ^ 2v/2a3/2 + (2ci + 0)^/2 - 

3{d — a) 

(2d + 0)3/2 + (2d - 2a)3/2 - (2d - 2a + 0)3/2 - (2a - 2d + 0)3/2). (23) 

The experiment for 70GeV protons in Si < 111 > gives the average reflection angle, 39.5 =F 
2.0/irarf Q. Using the data = = 0.58 * 10-^ a = 0.78A, d = 3.13A, the formula ([23]) 
gives x+ = 2 * a+ = 37.6firad. This result is lower than experimental, but it stays inside 
the interval of experimental accuracy. 

Now, we estimate the "refracted" part of the angle distribution for positive particles. The 
integral for average angle will be 



-^^ Hi-d)v^ 6(a/1-62_ yi_0_52) 

a_ = — 1= / ( 

aVl - J{i-d-d)VT^ V 1 - 



The estimation of this integral is also straightforward, but due to cumbersome result we 
write down only its infinitesimal approximation of the first order: 

«_ = 4(-203/2 _ 2^^3/2 ^ 4^2^/3/2 + (2a + 0)3/2 + 
3a 

(_2a + 0)3/2 - (2d + 2a)3/2 - (2d - 2a)3/2 - 
(2(i+0)3/2 + (2a + 2d + 0)3/2 + (-2a + 2rf + 0)3/2). (25) 

The same experimental parameters for 70GeV protons in Si < 111 > gives the negative 
average "refraction" angle X- = 2 * a;_ = —lO.Sfirad. We may guess here that the fraction 
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of particles on the right side of point B (see the emulsion's picture in [l|) may be treated as 
a volume "refraction" effect. It is separated from the channeling fraction and has a smaller 
average angle than the "reflection" angle. However, a proof (or disproof) of this is desirable. 

Now let us turn to the scattering of negative particles. The deflection function is described 
by the same equation (1191) with < 0. As a function of impact parameter b, it will cause the 
mirror reflection relatively the axis Ob and shift by the value ^/T^^^ in the direction of b 
(see Fig. 4(C), (D)). The corresponding fractions of scattering change its sign. In particular, 
the average "reflected" and "refracted" angles are given by the following equations: 

a+ = — (2(-0)3/2 + 2v^a3/2 - 4^2^3/2 _ _ ^.3/2 
3a 

(-2a - 0)=^/' + {2d + 2a)3/2 + (2d - 2df/^ + 
{2d - - (2a + 2d - (f)f^ - {-2d + 2d - 0)=^/'). (26) 

a. = -J- (-2(-0)3/2 + 2^2^3/2 _ 27253/2 _ (2d - 0)3/2 + 

3{d — d) 

(2d - 0)3/2 - (2d - 2d)3/2 + (2d - 2d - 0)3/2 + (2d - 2d - 0)3/2). (27) 

For negative particles, the "refraction" part is dominative and average angle of "refrac- 
tion" ( 1271) for 70GeV antiprotons in Si < 111 > is X- = 2 * a_ = — 37. 6/irad. The 
"reflective" part of the antiproton beam fl26l) will have small but positive average angle 
X+ = 2 * = lO.Sfirad. 

CONCLUSION 

We have found the deflection function (Eqs. f|T6|) . f|T7j) . f|T9l) ) for relativistic scattering by 
the system of periodic rectangular rings. Analysis of this function shows that an incident 
collimated beam is divided by each crystal's period into two fractions: "reflected" and 
"refracted" . These fractions are separated by a gap depended on potential shape, range of 
energy or radius of the bend. The particles in the "refracted" fraction are mostly deflected 
to the side of the ring's bend. We receive analytical expression for average "reflected" 
angles (Eqs. (1231) . (!26|) ) which seems in satisfactory agreement with experiment. We also 
receive expression for average "refracted" angles (Eqs. fl25|) . fl27|) ). The "refracted" fraction 
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is smaller than "reflected" for positive particles, but it might be present in scattering data. 
This fraction should dominate in the scattering of the negative particles. 
Author is grateful to Yuri M. Ivanov for the technical report j^. 
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